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A BAYESIAN  MODEL  FOR  DETERMINING  THE  OPTIMAL 
TEST  STRESS  FOR  A SINGLE  TEST  UNIT 


By 

H.  F.  Martz,  Jr,  and  M.  S.  Waterman 


ABSTRACT 

Consider  the  case  of  a single  test  unit  which  must  be  tested  at  some 
level  of  test  stress.  Suppose  that  the  test  stress  level  is  free  to  be 
determined,  and  that  only  the  survival  or  nonsurvival  of  the  unit  is 
observed.  It  is  assumed  that  the  unit  is  designed  to  withstand  a known 
and  specified  design  stress  level.  A Bayesian  model  is  developed  for 
determining  the  required  level  of  test  stress  which  maximizes  the  expected 
probability  of  survival  at  the  design  stress  .level . Engineering  experience 
from  similar  past  tests  on  similar  units  is  used  to  fit  the  model,  A 
practical  application  illustrates  the  method.  The  sensitivity  of  the 
procedure  to  changes  in  the  parameters  used  in  fitting  the  model  is  also 
examined.  The  procedure  is  fairly  insensitive  to  throe  parameters 
required  in  fitting  the  model  in  the  example.  . 

I.  INTRODUCTION  l\ 

In  many  engineering  testing  situations,  only'^a  single  unit  is  available 
for  testing.  The  unit  may  be  a component,  subsystem  or  complete  system. 
Suppose  that  a single  test  unit  is  to  be  tested  at  some  level  of  a single 
test  stress  which  is  free  to  be  selected  by  the  test  engineer.  Further, 
suppose  that  the  unit  is  designed  to  withstand  a known  and  specified  level 
of  design  stress.  It  is  further  assumed  that,  once  the  test  is  conducted, 
only  the  survival  or  nohsurvival  of  the  unit  is  observed.  i •• 
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For  example,  consider  the  case  of  a fuel  container  of  a radioisotope 
thermoelectric  generator  (RTG)  system,  which  is  the  power  supply  for  a 
space  satellite.  The  radioisotope  fuel  container  is  designed  to  withstand 
a certain  impact  onto  flat  plate  steel,  such  as  might  occur  during  a launch 
pad  overpressure  accident.  As  part  of  the  required  safety  analyses, 
tests  designed  to  simulate  such  an  accident  must  be  performed.  A proto- 
typic  unit,  using  simulated  fuel,  is  impacted  onto  flat  plate  steel  at 
some  test  velocity  which  must  be  determined.  This  test  velocity  may  or 
may  not  be  taken  to  be  the  design  velocity.  According  to  a precise  defini- 
tion of  failure,  e.g.,  if  the  unit  ruptures  to  the  extent  that  one  or  more 
fuel  elements  are  exposed,  the  unit  either  survives  or  fails  the  test. 

This  example  will  be  further  considered  in  Section  IV. 

The  model  to  be  developed  incorporates  the  following  aspects.  Suppose 
that  the  test  unit  survives  a test  stress  which  exceeds  the  design  stress. 
It  is  reasonable  that  this  should  increase  the  experimenter's  confidence 
in  the  ability  of  similar  units  to  survive  the  design  stress.  On  the 
other  hand,  if  the  test  unit  is  tested  at  too  great  a stress,  the  unit  is 
likely  to  fail,  thus  providing  little  information  about  the  unit's  ability 
to  survive  the  design  stress.  The  model  effectively  trades  between  these 
two  alternatives  in  seeking  the  optimum  desired  level  of  test  stress 
The  precise  definition  of  "optimum"  will  be  discussed  in  the  next  section 

The  philosophy  of  the  proposed  model  is  to  test  at  a high  enough 
stress  level  to  provide  assurance  but  not  failures.  This  contrasts  with 
the  usual  statistical  philosophy  which  is  to  test  at  various  levels,  some 
of  which  are  high  enough  to  insure  failures.  Of  course,  more  than  one  unit 
must  be  tested  in  this  case.  Easterling  (1975)  develops  an  over-test 
procedure,  referred  to  as  a "sensitivity  test",  which  is  based  on  such  a 
statistical  philosophy.  Much  of  the  literature  on  accelerated  life  testing 
considers  the  effect  of  stress  on  certain  failure  characteristics.  An 
excellent  bibliography  on  accelerated  life  testing  is  provided  by  Lowe  and 
Waller  (1975). 

II.  THE  MODEL 

Let  Sj^  denote  the  event  that  the  test  unit  survives  a test  of  stress 
k-S^,  where  is  the  given  design  stress.  Also  let  Pj^  = Prob(Sj^).  A 
Bayesian  approach  is  used,  in  which  the  uncertainty  in  Pj^  is  e.xprcsscd  by 
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assuming  that  is  a random  variable  having  a modified  negative- log  gamma 
prior  distribution  with  probability  density  function  (pdf)  given  by 


f (p;a,6,5) 


p(l/ekh-l 

BV‘r(a) 


0<P<1, 

0<k£“, 

0<a,B,6 


(1) 


This  distribution  may  be  derived  from  the  fact  that,  if  X has  a gamma 

distribution  with  shape  parameter  a and  scale  parameter  6,  then  T=exp[-Xk  ] 

has  the  distribution  given  in  (1).  Here  k is  the  test  stress,  expressed 
in  units  of  design  stress  S^.  The  parameter  6 appearing  in  (1)  is  used  to 
rescale  k,  for  reasons  to  be  discussed  in  the  next  section.  The  usual 
negative-log  gamma  distribution  may  be  obtained  by  letting  k=6=l.  The 
negative-log  gamma  distribution  has  been  previously  discussed  and  used  in 
reliability  by  Springer  and  Thompson  (1965,1967),  Mann  (1970),  and  Mastran 
and  Singpurwalla  (1974). 

The  mean  and  variance  of  (1)  are 

E(P;a,6,6)  = (1+Bk‘^)'“  (2) 

and 

V(P;a,B,6)  = (l+2Bk‘^)"“  - (l+6k'^)‘^“  , (3) 


respectively.  It  might  be  expected  that  the  mean  survival  probability 
curve  have  a reflected  S-shape  as  a function  of  k.  The  mean  given  in  (2) 
has  this  property  for  certain  combinations  of  a,B,6.  It  is  illustrated 
in  Figure  1 for  several  choices  of  a with  B and  6 computed  according  to 
the  example  in  Section  IV,  Figure  2 shows  the  standard  deviations  for 
the  same  set  of  distributions.  A procedure  for  identifying  a,B  and  6 
will  be  presented  in  the  next  section. 

The  distributions  of  interest  are  the  two  posterior  distributions 
of  P,  conditional  on  the  survival  (nonsurvival)  of  the  test  unit.  By  a 
simple  application  of  Bayes'  Theorem,  we  obtain  the  two  posterior  ]idf's 


f(p  |s^;a,B,6) 


B“k“‘^r(a)(l  + Bk'^)"“ 


0<p< 1 , 

a,B,5>0 


(4) 
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, Cl/Bk  )-l,  p ,a-l 
^ 3 k“*r(a) [l-(l+3k  ) “] 


0<p<l, 

0 < k <oo , 

a ,3 ,6>0 , 


(5) 


where  denotes  the  event  that  the  test  unit  does  not  survive  a 

test  of  stress  k-S  . The  cumulative  distribution  functions  (cdf’s) 

o 

associated  with  (4)  and  (5)  may  be  expressed  in  terms  of  the  chi- 
2 

square  (x  ) distribution  as 

F(p|Sj^;a,S,6)  = Prob  j P<^p  | Sj^;a,B  .sj  = Prob 

and 

F(piSj^;a,  3,6)  = Prob  | P<p  j ;a  , 3 , 6 | 


v2  . -2(1-H3k°)fnp 
6k* 


Prob 


2 -2£np 


- (l+Bk*^)  ’“Prob 


X?  . -2(1^3k’^)£np 


^2a  gj^6 


1-  d + 3k'^) 


(7) 


2 2 

where  denotes  a x random  variable  with  2a  degrees  of  freedom, 


The  posterior  means  of  (4)  and  (5)  are  easily  computed  to  be 

-a 

a + Bk' 


E(P^|S^;a,6,6) 


(8) 


and 


E(Pj^lS^;a,3,6)  = 


(l+Bk*^)""  - (l+2Bk'^)"“ 


l-Cl  + Bk*^)"" 


(9) 


The  unconditional  probability  of  survival  and  nonsurvival 
of  the  test  unit,  when  tested  at  stress  level  k,  are 


and 


Prob(Sj^)  = (1  + Bk^)'“  , 


Prob(S^)  = 1-(1+Bk^)'“  , 


(10) 


(11) 


respectively.  Note  that,  since  only  a single  unit  is  to  be  tested, 
(10)  is  the  same  as  (2). 


(6) 


figure 
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Thc  model  for  use  in  obtaininj’  the  required  test  stress  level 
is  based  on  the  following  two  propositions: 

. If  the  test  unit  survives  a test  stress  which  exceeds  the 
design  stress,  then  remaining  units  should  have  a higher 
expected  probability  of  survival  at  the  design  stress  than 
if  the  test  unit  had  been  tested  and  survived  at  the  design 
stress . 

. If  the  test  unit  does  not  survive  a test  stress  which  exceeds 
the  design  stress,  then  remaining  units  should  have  either 
t-he  same  or  a higher  expected  probability  of  survival  at  the 
design  stress  than  if  the  test  unit  had  been  tested  and  failed 
at  the  design  stress. 

In  both  propositions,  the  increase  depends  upon  the  difference 
between  the  test  and  design  stress  levels.  Mathematically,  let  us 
quantify  the  first  proposition  above  according  to 

E(PllSj^)  = g^(k)E(PjS^),k>l,  (12) 

where  g]^(k)  is  a suitably  chosen  function  of  k which  has  the  follow- 
ing properties:  (i)  g^(l)  = 1;  (ii)  gj^(k)->-l/E(P^  | ) as  k->«>;  and 

(iii)  g^(k)^0,  at  all  points  of  continuity  of  g^(k).  The  second 
property  guarantees  that,  if  the  test  unit  survives  a test  of  infinite 
stress,  then  the  remaining  units  are  expected  to  survive  the  design 
stress  with  probability  equal  to  1. 

Similarly,  the  second  proposition  may  be  quantified  as 

lUPlU\)  = g^Ckll-CP  jSj)  ,k>l  , (13) 

whore  g2(^)  suitably  chosen  function  of  k which  satisfies  tlie 

following  properties:  (i)  g2(^)  ~ g2  (k) -’•H  (Pj ) /E  (Pj  I as 

and  (iii)  g2(l<)iP.  1 1 points  of  continuity  of  g,(k).  The 
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second  property  insures  that,  if  the  test  unit  fails  a test  of 
infinite  stress,  then  nothing  additional  has  been  learned  from  the 
test  about  the  ability  of  remaining  units  to  survive  the  design  stress 
Particular  choices  for  gj^(k)  and  g2(^^)  will  be  considered  in  Secs.  Ill 
and  IV. 

Now  let  us  consider  the  optimization  model  itself.  We  wish  to 
determine  the  value  of  k which  maximizes  the  expected  probability  that 
occurs,  given  that  we  test  a single  unit  at  test  stress  k.s^.  That 
is,  we  wish  to  maximize 


E(P^;a,6,5)  = E (P^  | ;a  , B , d)  Prob  (Sj^) 

+E(PilSk;a,B  d)Prob(Sj^)  . (14) 

Upon  substituting  (8) -(13)  into  this  expression  and  simplifying, 
we  obtain 

E(P^;a,B,5)  = [ Y^gi  (k)  - 7382  + Y 7^2  ’ 

where 

’1 

and 


(15) 

(16) 


\l+26/ 


a 


1-Y 


a 


1_ 

a 


(1+B)  - 1 (1+8)  - 1 

In  order  to  maximize  E(Ppa,B,5),  it  is  useful  to  solve 


3E(P^;  a,  6.  6) 
9k 


[Yigi'(k)-Y2?2^k)](l+Bk'^)  “ 


(17) 


,6-  1 

-a86[Y,g,(k)  - Yogot*^^!  L + Yo8^(k)  =0.  (18) 

(l+6k°) 

The  solution  to  (18)  yields  the  desired  optimal  test  stress  k^. 
The  solution  to  the  problem  of  finding  k such  that  (14)  ir  maximized 
is  discussed  in  the  next  section. 
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III.  FITTING  THE  MODEL 

First,  let  us  consider  a procedure  for  estimating  a,  3,  and  6.  Since  only 
a single  test  unit  is  assumed  to  be  available,  the  suggested  procedure  is  neces- 
sarily subjective.  Consider  the  following  two  questions; 

Question  1:  Prior  to  the  test,  at  what  stress  level  will  the  test  unit 

have  approximately  a 95%  expected  chance  of  survival? 

Question  2:  Prior  to  the  test,  at  what  stress  level  will  the  test  unit 

have  approximately  a 5%  expected  chance  of  survival? 

Now  k2>k|  and  both  are  expressed  in  units  of  design  stress.  Then 


(1+B  = C.,  i=l,2, 


where  = 


.95  and 


.05.  Then 


3 = i=l,2. 


and  further  simple  algebra  yields 


5 = 


ln[(C’^/“-l)/a‘^/°‘-l)] 
rni-kY/kp"  


(18) 


(19) 


Therefore  we  have  3=3 (a)  and  6=6 (a),  so  that  our  three  parameter  model  has  been 
reduced  to  a one  parameter  model . 

Since  the  prior  variance  of  the  survival  probability  is  given  by 

VCp;a,6,6),  = V(a)  = (l+23k‘^)‘“-  Cl+Bk'^)"^“, 

the  parameter  a may  be  chosen  to  coincide  with  the  experimenter's  prior  estimate 
of  the  variation  at  some  stress  k.  For  the  examples  we  have  worked  out,  V(a)  has 
been  observed  to  be  a decreasing  function  of  a. 

Two  other  functions,  g^^(k)  and  g2(k),  must  be  specified.  In  our  calcula- 
tions, we  have  taken 

g^(k)  = k"",  k<K 

and 

g2(k)  = 1,  k<K 

where  K is  an  unspecified  upper  limit  of  test  stress  beyond  the  range  of  practical 
interest.  For  k>K,  suitable  adjustments  would  have  to  be  made  to  these  choices 
of  gj(k)  and  g^(k)  to  ensure  that  the  appropriate  asymptotic  properties  discussed 
in  .Section  IT  arc  present. 
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Recall  that 


ECPjtSk)  = gjCk)  ECPjsp 


= E(Pjsp 

The  choice  g2(k)  = 1 expresses  the  situation  in  which,  if  the  test  unit 

fails  a test  of  stress  k,l<k<K,  then  the  expected  probability  of  survival  is  the 

same  as  if  the  test  unit  had  failed  a test  at  the  design  stress. 

To  determine  c,  consider  the  following  question: 

Question  3:  Prior  to  the  test,  suppose  that  a hypothetical  test  unit  was 


tested  and  survived  an  increased  test  stress  level. 

ic 


At  what  stress  level  k^  will 


the  expected  lailu/s  probability,  E(l-P.,  !s.  ),  be  one-half  as  large  as  the  expected 
failure  probability  of  a hypothetical  unit  which  was  tested  and  survived  the  design 
stress,  E(l-Pj|s^)? 

Of  course,  must  also  be  expressed  in  units  of  design  stress.  Then 

1 - E(P|Sj)  = 2 I 1 - E(PlSk  )| 


or 

Some  elementary  algebra  yields 


-a 


= 2 


*^3  \l+6  / 


In 


/ln[k3] 


C20) 


Therefore,  for  the  specific  choice  of  gj(k)  and  g2(k),  the  parameters  6,6,  and 
c are  given  by  equations  (18),  (19),  and  (20).  Then,  when  a is  chosen  in  accord 
with  the  experimenter's  estimate  of  the  variation,  the  parameters  of  the  model  are 
completely  determined. 


IV.  EXAMPLE 

As  indicated  in  the  introduction,  the  example  considered  here  concerns  the 
fuel  container  of  a radio-isotope  thermoelectric  generator  (RTG)  system  used  as 
the  power  supply  for  a space  satellite.  In  the  radioisotope  fuel  container  there 
are  a number  of  fuel  elements,  which  are  simply  spheres  which  contain  the  fuel. 
The  RTG  is  designed  to  withstand  impact  onto  flat  plate  steel  at  a certain  veloc- 
ity so  that,  for  example,  launch  pad  accidents  will  not  release  radioactive 
material  to  the  environment. 
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The  answers  to  questions  1,  2,  and  3 of  Section  III  were  solicited  from  a 
group  of  engineers  at  the  Los  Alamos  Scientific  Laboratory  responsible  for  such 
impact  tests.  The  answers  for  one  particular  RTG  S’/stem  of  interest  were  as 
follows: 

1)  At  approximately  140  fps,  the  test  unit  should  have  roughly  a 95% 
expected  chance  of  survival. 

2)  At  approximately  180  fps,  the  test  unit  should  have  roughly  a 5% 
expected  chance  of  survival. 

5)  At  approximately  135  fps,  E(l-Pj|Sj^  ) should  be  roughly  one-half  of 
E(l-PjSj).  ^ 

The  design  stress  velocity  is  100  fps.  Thus,  kj=1.4,  k2=l-8,  and  k2=1.35. 
Equation  (14),  for  g^  and  g2  specified  in  section  III,  becomes 

e(p^;  a, 3 (a), 5 (a))  = (l+3k*^)'°‘ 

/ \ / \ 

The  optimal  k,  k^,  is  found  on  a computer  by  a simple  search  program. 

Certain  a and  k ,k-,  k_  result  in  difficulties  in  computation. 

For  example,  if  8k°=0,  then  on  the  computer 

1 - (l+6k’^)"“  = 0 

whereas  a simple  expansion  shows  that 
1 - (l+Bk®)"^^  = + aek*^ 

IVhile  such  approximations  were  used  whenever  possible,  it  was  not  possible  to 
obtain  all  values  in  the  tables  below.  A indicates  that  the  computation  was 
not  performed. 

Table  1 gives  the  values  of  3=3(a),  6=6(a),  c=c(a),  and  k^=k  (a)  for  several 
different  choices  of  a.  The  values  of  are  kj^=.14,  k2=1.8,  and  k^^l.35. 

Observe  that  k^  is  a very  stable  function  of  a. 
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TABLE  1 

VALUES  OF  6,  6,  c,  AND  OPTIMAL  STRESS  k FOR  SELECTED  VALUES  OF  a 

o 


a 

B 

6 

c 

k 

o 

.3 

3.07xl0‘® 

46.42 

1.53x10"® 

1.96 

.4 

4.22x10"^ 

37.71 

2.81x10"^ 

1.35 

.5 

.0000018 

32.69 

.0000015 

1.36 

.6 

.0000044 

^ 29.45 

.0000044 

1.37 

.7 

.0000080 

27.23 

.0000093 

1.38 

.8 

.0000120 

25.61 

.0000159 

1.38 

.9 

.0000160 

24.38 

.0000240 

1.39 

1.0 

.0000198 

23.43 

.0000330 

1.39 

2.0 

.0000370 

19.48 

.0000123 

1.43 

3.0 

.0000365 

18.30 

.0001825 

1.45 

4.0 

.0000330 

17.74 

.0002198 

1.46 

5.0 

.0000294 

17.42 

.0002450 

1.47 

10.0 

.0000181 

16.79 

.0003021 

1.48 

15.0 

.0000129 

16.58 

.0003232 

1.49 

20.0 

.0000100 

sl6.48 

.0003342 

1.49 

25.0 

.0000082 

16.42 

.0003409 

1.49 

30.0 

.0000069 

16.38 

.0003454 

1.49 

35.0 

.0000060 

16.35 

.0003487 

1.49 

40.0 

.0000053 

16.33 

.0003512 

1.49 

45.0 

.0000047 

16.32 

.0003530 

1.49 

50.0 

.0000042 

16.30 

.0003546 

1.49 

Values  of  the  prior  standard  deviation. 

1/2 

V ' (a) , are  plotted 

as  a function 

of  stress  k 

in  Figure  2 for  several  choices  of  a.  As  can  be  seen. 

the  prior 

standard  deviation  decreases  as  a 

increases . 

Figure  3 gives  a plot  of  E(P), 

E (P  ]S^)  and 

E(P|Sk)  as  a function 

of  stress  k 

for  a=0.1.  It  is  observed  that 

i'Cn|s)  < Ecn  < F.(pjsi . 


FTQ(RF. 


FinURE 


FIGURE  7 
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Also,  it  is  easy  to  see  that  the  asymptote  for  E(P|S)  is 


lim 

k-«o 


E(P|S)  = lim 
k-x» 


il2^ 

1+Bk® 


2 


-a 


Figures  4-7  give  plots  similar  to  Figure  3 for  a = 0.5,  1.0,  10.0.  and 
50.0,  respectively. 

Earlier,  in  section  III,  we  suggested  the  a could  be  chosen  to  coincide 
with  the  experimenter's  prior  estimate  of  the  variance  at  a given  stress  k. 

It  is  also  possible  to  choose  a by  use  of  Figures  3-7.  From  these  figures, 
it  is  observed  that  the  difference  between  the  prior  and  posterior  expected 
survival  probabilities  is  larger  for  smaller  values  of  a.  That  is,  for  small 
values  of  a,  the  expected  survival  probability  is  more  sensitive  to  the  test 
result  than  for  large  values  of  a. 

Let  us  now  examine  the  sensitivity  of  the  optimal  test  stress  k^,  as  given 
by  the  solution  to  (21),  to  variations  in  the  answers  to  questions  1,  2,  and 
3.  It  is  important  to  do  this  since  the  answers  to  these  questions  may  be  in- 
accurate. Such  inaccuracy  may  be  due  to  either  lack  of  precise  knowledge  by 
the  person(s)  answering  the  questions  or  lack  of  clear  understanding  of  the 
precise  information  being  solicited  in  the  questions. 

First,  consider  the  sensitivity  of  the  optimal  test  stress  to  changes 
in  ki  and/or  kj.  Tables  2-5  give  the  optimal  test  stress  ko  as  a function  of 
several  choices  of  kj  and  k^  for  a = 0.5,  1.0,  10.0,  and  50.0,  respectively. 

In  Tables  2-5,  k^  = 1.35. 


TABLE  2. 

OPTIMAL  TEST  STRESS  k FOR  SEVERAL 
VALUES  OF  kj  AND  k^  (a  =°0.5  AND  k^  = 1.35) 


1.2 

1.3 

1.4 

1.5 

1.6 

1.6 

1.17. 

1.26 

1.87 

* 

•k 

1.7 

1.18 

1.27 

1.36 

* 

* 

1 .S 

1.18 

1.27 

1 . 36 

★ 

k 

1 

1.1!’ 

1.27 

1 ..57 

•k 

1.90 
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TABLE  3. 
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TABLE  5. 

OPTIMAL  TEST  STRESS  kg  FOR  SEVERAL 
VALUES  OF  kj  AND  k^  (a  = 50.0  AND  k^  = 1.35) 


It  is  observed  that  the  optimal  test  stress  ranges  between  1.09  (a  = 50, 
kj  = 1.6,  kj  = 1.7)  and  1.93  (a  = 0.5,  k^  = 1.6,  k^  = 2.0).  For  a given  a, 
the  optimal  test  stress  is  fairly  insensitive  to  changes  in  for  small 
values  of  k^.  On  the  other  hand,  for  a given  a,  the  optimal  test  stress  is 
fairly  insensitive  to  changes  in  kj  for  large  values  of  k^ . As  both  kj  and 
kj  increase,  the  optimal  test  stress  is  fairly  stable. 

Now  consider  the  sensitivity  of  the  optimal  test  plan  to  changes  in  k^, 
since  this  quantity  was  held  fixed  in  Tables  2-5.  Table  6 gives  the  optimal 
test  stress  k^  as  a function  of  several  choices  of  k^  and  a for  the  nominal 
values  kj  = 1.4  and  k^  = 1.8.  It  is  observed  that  k^  is  quite  insensitive  to 


TABLE  6. 

OPTIMAL  TEST  STRESS  Lq  FOR  SEVERAL 
VALUES  OF  k3  AND  a (k,  = 1.4  AND  k^  = 1.8) 


1 

j 

1.10 

1.20 

1.30 

1.35 

1 .40 

1 .50 

0.5 

1.39 

1..38 

1 . 36 

1.36 

1.36 

1 . 35 

1.0 

1.43 

1.41 

1.40 

1.39 

1.39 

1 . 38 

a 

10.0 

1.49 

1.49 

1.48 

1.48 

1.48 

1 .48 

.50.0 

1.50 

1.50 

1.49 

1.49 

1 .49 

1 .49 
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changes  in  kj  for  a fixed  value  of  a.  This  is  an  important  result,  since  the 
answer  to  Question  3 is  likely  to  be  somewhat  arbitrary  in  practice.  That 
is,  in  practice,  kj  may  be  an  imprecisely  known  value. 

V.  CONCLUSIONS 

A Bayesian  procedure  for  determining  the  optimal  test  stress  for  a single 
test  unit  has  been  developed.  The  procedure  is  both  objective  as  well  as  sub- 
jective. It  is  subjective  in  the  sense  that  the  necessary  parameters  in  the 
model  are  estimated  from  best  available  information  prior  to  the  test  results. 
These  estimates  are  then  used  in  an  objective  manner  to  provide  the  required 
test  stress.  The  test  stress  provided  by  this  procedure  is  "optimal"  within 
the  model  framework  in  a certain  well-defined  sense.  Namely,  this  optimal 
test  stress  maximizes  the  modeled  expected  unconditional  probability  of  survival 
at  the  design  stress.  The  model  effectively  trades  between  two  extremes.  The 
first  represents  the  increasing  likelihood  of  survival  at  the  design  stress 
gained  as  a result  of  a test  unit  surviving  increasing  test  stress.  This  gain 
is  countered  by  the  correspondingly  decreasing  probability  of  test  unit  sur- 
vival as  the  test  stress  increases. 

The  model  was  used  to  determine  the  impact  test  velocity  in  an  impact 
test  of  a certain  radioisotope  fuel  container.  The  optimal  test  velocity  was 
found  to  be  approximately  30-50  percent  above  the  design  impact  velocity.  In 
addition,  a limited  sensitivity  analysis  to  the  subjective  estimates  required 
in  fitting  the  model  was  conducted.  It  was  observed  in  this  example  that  the 
optimal  test  velocity  was  fairly  insensitive  to  the  subjective  estimates. 

This  may  or  may  not  be  true  in  other  applications.  Consequently,  as  a safe- 
guard, it  is  recommended  that  such  a sensitivity  analysis  be  routinely  con- 
ducted when  applying  this  model  in  practice. 
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